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We study the energy spectrum of symmetric double quantum dots in narrow-gap carbon nanotubes 
with one and two electrostatically confined electrons in the presence of spin-orbit and Coulomb 
interactions. Compared to GaAs quantum dots, the spectrum exhibits a much richer structure 
because of the spin-orbit interaction that couples the electron's isospin to its real spin through two 
independent coupling constants. In a single dot, both constants combine to split the spectrum 
into two Kramers doublets, while the antisymmetric constant solely controls the difference in the 
tunneling rates of the Kramers doublets between the dots. For the two-electron regime, the detailed 
structure of the spin-orbit split energy spectrum is investigated as a function of detuning between the 
quantum dots in a 22-dimensional Hilbert space within the framework of a single-longitudinal-mode 
model. We find a competing effect of the tunneling and Coulomb interaction. The former favors a 
left-right symmetric two-particle ground state, while in the regime where the Coulomb interaction 
dominates over tunneling, a left-right antisymmetric ground state is found. As a result, ground 
states on both sides of the (ll)-(02) degeneracy point may possess opposite left-right symmetry, 
and the electron dynamics when tuning the system from one side of the (ll)-(02) degeneracy point 
to the other is controlled by three selection rules (in spin, isospin, and left-right symmetry). We 
discuss implications for the spin-dephasing and Pauli blockade experiments. 

PACS numbers: 73.63.Fg,71.70.Ej,73.21.La,73.63.Kv 



I. INTRODUCTION 

Coherent control over the charge and/or the spin of an 
electron or hole is a key ingredient for quantum compu- 
tation or spintronic devices. It is of importance to have 
coupled two level systems (qubits) that can be controlled 
and manipulated efficiently without loss of the stored in- 
formation. A promising candidate and natural two-state 
system for a robust qubit is the spin of an electron. Co- 
herent manipulation as well as preparation and read-out 
of a single confined spin in few electron semiconductor 
quantum dot (QD) systems have been demonstrated, see 
Refs. [T] and [5] and references therein. 

Spin qubits in carbon nanotubes^' (CNT) are believed 
to be even more robust due to the absence of hyper- 
fine coupling in 12 CP This is in contrast to GaAs QDs 
where the phase coherence suffers from hypcrfine cou- 
pling due to the nuclei of the host crystal. However, 
CNTs pose other challenges and complications. First, 
few electron QDs are not easily fabricated and, secondly, 
the isospin degree of freedom present in the honeycomb 
carbon lattice provides another quantum two level sys- 
tem that must be included in the analysis. 

The band structure of electrons in nanotubes can be 
understood starting from that of graphene which has a 
linear dispersion relation similar to massless Dirac-Weyl 
fermions. Graphene is a zero gap semiconductor, but 
when the graphene sheet is rolled to form a nanotube, the 
quantization condition for nanotubes leads to metallic or 
semiconducting behavior, depending on chiralityP^ The 
curved geometry creates a mass term in the Dirac spec- 



trum and thus a bandgap even for the nominally metallic 
tubesP^HHI This bandgap allows for electrostatic confine- 
ment of electrons and creation of few electron QDs, other- 
wise not possible due to the Klein paradox.^ Recent ex- 
periments have shown that it is indeed possible to confine 
electrons in single^US] anc j double quantum dots (DQD) 
in a CNT by means of ele ctrostat ic gates in cleanly grown 
small bandgap nanotubes! 5 * 15 * 20 ! The present study is mo- 
tivated by these experimental results. 

In the recently investigated few-electron nanotube 
QDs, the four-fold degeneracy due to the spin and isospin 
degrees of freedom is split by SO coupling, giving rise to a 
coupling of spin and isospin degrees of freedom. In plane 
graphene, spin-orbit (SO) coupling (being a relativistic 
effect) of its 7r-electrons is of dozens of /xeV onljl^HUl 
and therefore of minor importance. The curved geometry 
of nanotubes induces SO coupling on the order of 10 _1 
meV among the single particle levels of the electrons.^ 
While the curvature- induced SO co upling was envisioned 
previously for semiconductors^ 5 -^ f or nanotub es it is the 
dominant mechanism. It was Ando^and others^ ^ who 
developed the first theories of SO coupling in nanotubes. 
More recent theoretical investigations extended this work 
by including the a- and 7T- bands in full as well as the 
curved bonds between neighboring atomsP^HUl Lowest 
order perturbation theory shows that the SO coupling 
is inversely proportional to the radius of curvature and 
originates from the intra-atomic SO coupling in a car- 
bon atom. Even though this is a weak coupling com- 
pared to heavier atoms, the combined effect of curvature 
and intra-atomic SO coupling splits a four-fold degen- 
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erate level into two Kramers doublets by a fraction of 
a meV. The demonstration of electrostatically con fined 
particles in CNT-QDs in the presence of SO coupling^ 18 * 20 ! 
has motivated recent theoretical investigations. The sin- 
gle electron QD setup and in particular the influence of 
the electron-phonon coupling on the decoherence are sub- 
jects of a work by Bulaev et. al.^ The two-particle prob- 
lem has been studied numerically in Refs. 1551 and l3"ol for 
a hard wall and a harmonic potential, respectively Fur- 
thermore, hyperfine interactions and t heir c onsequences 
for Pauli blockade have been discussed! 37 * 38 * 

Here, we present a theory of the energy spectrum for 
a CNT-DQD in the presence of SO coupling in the enve- 
lope function formalism within an exact diagonalization 
scheme. The system is described by three quantum num- 
bers [spin s, isospin t, and left and right dots (L/R)], 
each taking two values, and by the discrete and con- 
tinuum spectra of the longitudinal motion. Depending 
on the specific parameter values of the DQD that vary 
in a wide range, two electrons confined in a DQD find 
themselves in rather different regimes. In this article, 
we concentrate on coated narrow-gap nanotubes of the 
kind investigated in Refs. l5l and [20l that seem well suited 
for spintronic applications. In such DQDs, a small elec- 
tron mass to* increases the separation between the levels 
of longitudinal quantization, while a high dielectric con- 
stant n ~ 10 suppresses the Coulomb repulsion, which 
would otherwise result in level mixing. Restricting our- 
selves with the lowest longitudinal mode, we concentrate 
instead on the detailed structure of the energy spectrum 
emerging from the spin-isospin coupling, and the depen- 
dence of the interdot tunneling and Coulomb energies on 
the symmetry of wave functions and SO coupling. How- 
ever, we stop short of discussing the influence of phonon 
and hyperfine couplings as well as the scattering between 
different isospin statesPS It should be mentioned that 
similar physics occurs in silicon double quantum dots 
where, however, the valley degeneracy is usually broken 
and spin-orbit interaction does not play an important 
rolePl 

To set up our model calculation, we use the well- 
established model for p 2 -orbitals of graphene to describe 
electrons in a nanotube. The DQD confinement is mod- 
eled by a double square-well potential along the axial 
direction of the nanotube. We take into account SO cou- 
pling effects on the single particle levels and discuss their 
influence on the energy spectrum in the presence of an 
axial magnetic field B. In the framework of a single- 
mode approximation, and using the eight-function basis 
(8 = 2 x 2 x 2) of single-electron states, we present a sym- 
metry classification of the eigenstates of a two-electron 
symmetric DQD in its 22-dimensional Hilbert space. We 
find its energy spectrum numerically for the comparable 
values of the tunneling integral, the Coulomb interaction 
and the SO splitting. Our main results include the de- 
tuning dependence of the spectrum and the effect of a 
magnetic field lifting all spectrum degeneracies. Finally, 
we discuss challenges and opportunities for experimental 



studies. 

The structure of the article is as follows. After we have 
summarized the physical properties of a single electron 
in a nanotube in Sec. |II A| we turn in Sec. |II B| to the 
model of an electrostatically generated DQD and solve 
the eigenvalue problem for a DQD with a square-well 
potential. We clarify the effect of two SO coupling con- 
stants, A and A l5 on the energy splitting between the 
Kramers doublets and on the tunneling integral. Sec. HI 
presents the symmetry classification of the two-electron 
wave functions as well as a description of the techniques 
used for calculating Coulomb integrals on SO modified 
wave functions, and the effect of SO coupling on Coulomb 
integrals as well as their Independence 
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In Sec 

present our main results on the energy spectrum of a two- 
electron DQD as a function of detuning and its transfor- 
mation when a magnetic field is applied. Summary and 
discussion are presented in Sec. |V| 



II. MODEL 

A. Single electron in a nanotube 

We consider a single-wall CNT whose electronic prop- 
erties are described within a tight-binding model for the 
p z -orbitals of neighboring carbon atoma^. As usually, 
we solve for the band structure of a plane graphene 
sheet first and then impose periodic boundary conditions 
for the electronic motion along the circumferential di- 
rection defined by a chiral vector C. It is defined as 
C = 77,iai + 77,2a2 in terms of the primitive lattice vec- 
tors ai = a(l,0) and &2 — a(l/2, -\/3/2), where (711,712) 
are integers. Coordinates along the circumferential and 
translational directions, C and T in Fig. [I] are (c, t) . Due 
to the honeycomb lattice structure of graphene, nearest- 
neighbor atoms belong always to different sub-lattices 
A and B; the lattice constant is a — 0.246 nm. The 
chiral angle, which is the angle between C and ai, is 
9 = arctan[v / 3772/(2ni + 772)]. In graphene, two spin- 
degenerate 7r-bands (the conduction and valence bands) 
cross at six vertices of the Brillouin zone. Two pairs 
of translationally nonequivalent vertices, K r , form two 
Dirac points; therefore graphene is a semi-metal. Here 
we choose the Dirac points as K T = 27r/a(r/3, 1/V3), 
with r = ±1. The effective low-energy Hamiltonian is 
obtained by expand ing in the electron momentum near 
the Dirac points K T pI™ 



% = Hv(T 3 k c ai + k t a 2 ) 



U) 



Here, T3 is the diagonal Pauli matrix with eigenvalues 
t = ±1 in the K r -isospin subspace. The Pauli matri- 
ces <jj (j = 1,2) act in sublattice space and account for 
the two carbon atoms in the primitive unit cell of the 
honeycomb lattice (up to a unitary transformation that 
depends on 9, see Ref . I34p . The quasimomentum compo- 
nents along the C and T axes are k c and k t , see Fig. [T] 
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FIG. 1: Sketch of the planar graphene sheet showing the 
honeycomb lattice structure. Nearest neighbors belong to two 
different sublattices A and B. A nanotube with chirality (4, 1) 
is formed when the sheet is rolled up along the direction of 
the chiral vector C. The chiral angle 9 gives the misalignment 
between the chiral vector and the primitive lattice vector ai. 
The direction perpendicular to the chiral vector defines the 
tube axis and is denoted T. Within the tight-binding ap- 
proximation, we assume the z-direction to be perpendicular 
to T and C, i.e., the p z orbitals stick out of the plane of the 
figure. The Brillouin zone of the honeycomb lattice is shown 
in the upper right corner. The Dirac points are denoted by 
K T , and we choose in our calculations two inequivalent points 
K T = 27r/a(r/3, l/v/3), which are marked in red. A graphene 
sheet folded into a nanotube is shown in the lower right cor- 
ner. 



The eigenvalues of the Hamiltonian H are readily ob- 
tained 



E kc . kt = ±hv^kl + kl (2) 

where i-solutions in Eq. ([2]) correspond to the conduc- 
tion and valence band, respectively, that are degenerate 
in the isospin quantum number r = ±1, and v ~ 8 x 10 5 
m/s is the Fermi velocity in graphene. Plane-wave type 
eigenf unctions for the Hamiltonian % are 



jK t t 



exp {i(k c c + k t t)} ( k oM,c/v 



(3) 

where the coefficients zT , , for the conduction or va- 

k c ,kt ,c/v 

lence band (denoted by the subscripts c and v) are 



rk r — ikt 



rk r — ikt 



VWTkl 



(4) 



The position vector for the electron is r = r(x, y) = 
(ccos# — isin(9,csin# + tcosO) and the radius of the 
tube is R = \C\/2tt — a^n\ + n\ + n\n 2 l2-K. For 



(77.1,77.2) = (n,0)/(n,n) a CNT is called zigzag/armchair 
like, other nanotubes are called chiral!^ Obviously, for 
zigzag tubes 9 = and {x,y) = (c, t). In what follows, 
we restrict ourselves with conduction band electrons and 
designate amplitudes as zj,^ k . 

Wavcfunctions of nanotubes are periodic in the circum- 
ferential direction, i.e., \P(r) = \&(r + C). Consequently, 
the wavenumber k c of electrons/holes in this direction is 
quantized by the condition (k + K r ) • C = 2nm with m 
being integer numbers and k = (k c ,k t ). Depending on 
the direction and magnitude of C, there are two types 
of solutions obeying the periodicity condition around the 
circumference. If we use 



c = cxp {(27r?r/3)(ni — n 2 )} 



(5) 



as well as (ni — n 2 ) = 3M + v with v — 0, ±1 and M 
integer, we obtain the quantization of the wavenumber 
around the circumference as k c — > k m = (m — vt/3)/R. 
Note that v = is always fulfilled for armchair tubes, 
but for zigzag tubes only if n\ = 2>M . For v = ±1, the 
envelope wavefunctions accumulate phase factors. 

In graphene, a classification of quantum states by K T 
is protected by the conservation of the crystal momentum 
k. In nanotubes, however, it is at first sight not so clear 
that isospin is a good quantum number, since for some 
metallic tubes the folding of the graphene bandstructure 
onto the first Brillouin zone of the translational nanotube 
unit cell, results in Fermi points at k t = for both isospin 
values (sometimes classified as "zigzag-like" nanotubes^ 
or in Ref.lHl chapter 4 as metal-1). Nevertheless, for these 
nanotubes the isospin quantum number is protected by a 
screw axis of the order iVs(ni, n 2 ) defined by a Diophan- 
tine equation, see Ref. [321 For the k t = point, screw 
rotations are equivalent to spatial rotations, hence, it fol- 
lows from Eq. ^ that such rotations produce phase fac- 
tors exp[(27rir/ci7Vg)(ni — 712)] having complex conjugate 
values for r = ±1. Therefore, r = ± states belong to 
complex conjugate representations. For the other class 
of metallic tubes ( "armchair- like" or metal-2) where the 
Fermi points are different and at ±27r/3T, isospin is pro- 
tected by momentum conservation. This should clarify 
the meaning of the isospin quantum number r for specific 
types of nanotubes. 

If we insert the allowed quantized k m values into the 
dispersion relation of Eq. ([2]), we see that there is a gap 
between the conduction and valence bands given by 2E g 
with Eg = Hv\m — vt/3\/R. For m = and v = ±1, this 
gap between the conduction and valence bands is about 
2E g w 360 meV for R = 1 nm. Such nanotubes are thus 
semiconducting, whereas the v = tubes are nominally 
metallic. The curvature, however, opens a small gapf^Hlll 
likely causing the measured gaps of order 10 — 50 meV 
in Refs. [TH] and [501 The curvature effects appear in the 
Dirac Hamiltonian as a mass term, and Eq. (!]) reduces to 
a one-dimensional effective Hamiltonian with k c modified 
as (for the lowest energy mode m = 0) 



rk g , kg 
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3R 



(6) 
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where the last term scales with tube radius 
1/-R 2 and the induced g ap as E curv — Hvk curv oc 
1 /i? 2 lioiil] whjig specific expressions for the gap are 
model dependent, the o rder of magnitude estimate for 
R of a few nanometers i j 10 l 14 l 31 l 



E CUIV {0) ~ (hva/R 2 ) cos W ~ 10 meV. 



(7) 



We also include a magnetic field B which points in the 
tube axis direction T and induces an Aharonov-Bohm 
flux <&ab = BirR 2 through the cross section of the tube. 
This further modifies the circumferential momentum as 
rk g -» rk g + $ab/(R®o), with $ = hc/e being the flux 
quantum. Therefore, the nonrelativistic circumferential 
momentum fc" r equals 



(8) 



Besides the orbital effect, the magnetic field also leads 
to a Zeeman term given by StgfisB, where St — ±1/2 
is the spin projection along the CNT axis, fiB is the 
Bohr magneton, and g ~ 2 is the bare electronic g- fact or. 
This yields an energy difference between the different spin 
species of the electron. In this paper, we only consider 
tubes with finite gaps allowing for electrostatic confine- 
ment of electrons, and pay special attention to the tubes 
with the curvature-induced gaps (narrow-gap nanotubes, 
v = 0). We therefore write the (electron/hole) dispersion 
relation as 



kf + S t gfi B B. 



(9) 



Now we introduce the SO coupling which was shown 
to be an importan t effe ct in the recent experiments on 
few electron QDs.LUU^ol j n g enera i ; the coupling of the 
electron spin to its orbital motion is a relativistic effect 
for electrons moving in external electric fields. Asym- 
metric confinements in semiconductor QDs (extrinsic SO 
coupling, see Ref. S3] and references therein) can also 
provide such a coupling which, however, is one order of 
magnitude smaller than the SO coupling constants re- 
ported in CNTs. The pioneering theoretical work on 
the curvature induced SO derives the low-energy effec- 
tive Hamiltonian from a first-order perturbation theory 
in the atomic SO coupling as well as in the curvatureP-2 
Due to curvature, there are nonzero overlaps between 
the p z and p x<y orbitals of neighboring atoms. Combined 
with the atomic SO coupling which produces transition 
matrix elements between different quantum states on the 
same atom, a spin depende nt cou pling between the ad- 
jacent A and B atoms arisesp32Hl More recent work^SSl 
has extended this approach and added to the low-energy 
effective Hamiltonian of 7r-electrons in CNTs a term that 
is diagonal in sublattice (A, B) space. The generalized 
SO Hamiltonian near the Dirac points is 



Hso = A1T3CT1S3 + A0S3T3, 



(10) 



with S3 being a diagonal Pauli matrix acting in spin 
space. According to Eq. (10), the electron spin is still 



a good quantum number for the single particle problem, 
and in what follows the eigenvalues of S3 are denoted as 
s = ±1. The two coupling constants A and Ai depend 
on the type of the tube. Both are inversely proportional 
to the radius R, and Ao depends on the chiral angle as 
Ao oc cos (30), while Ai does not depend on 9. The term 
in Eq. |l0| ) proportional to Aq is diagonal in sublattice 
space and is responsible for the difference in the electron 
and hole spin-orbit gaps observed experimentally)^ see 



Eq. (17) below. The combination of the SO interaction 



with the effective nanotube Hamiltonian described above 
gives the final expression for the dispersion relation in the 
conduction and valence bands 



kl + k 2 + I A r + \ gi i B B ) s, (11) 



with 



9 R$o 



sAi 
hv 



(12) 



where the last term is a SO correction to k™ of Eq. (JsJ) . 
We note that k c is spin and isospin dependent, k c = k C:T>s , 
but to simplify notations we suppress the indices (r, s) in 
what follows. 

We note that the curvature and SO corrections to k c 
should be applied when calculating energy levels and 
Zk c ,k t amplitudes of Eq. (|3j). However, the phase factors 
of the wavefunctions ^ fe are fixed strictly by the pe- 
riodicity condition and cannot be changed by the renor- 
malization of k c nor by minor changes in R and K T due 
to the deformation of graphene when folding it into a 
nanotube. Finally, phase factors of the circumferential 
wavefunctions can be chosen quite generally as 



exp 



tM + n 2 ~- 



1 



(13) 



for the lowest circumferencial modes, m — 0, where tp 
is the azimuth along the circumference C, and r c is the 
circumferential component of r. In Eq. (13), the two 



first terms in the parentheses stem from the dot products 
(K r • C) and the v term is canceled by the leading term 
in k c of Eq. (|6). This equation will be used in Sec. MB 
below to calculate Coulomb integrals. 

According to Eq. ( |TT| ), the spin degeneracy of the 
single-particle levels is lifted by SO coupling that in the 
absence of a magnetic field opens a gap between the spin 
s =\\) and s =\\) states of size A.go(kt) for each of the 
K r points. It is defined as 



E 



- E 



+1 

k. c ,kt 



E 



-I 

fccife 



- E 



-t 

k c ,k t 



(14) 



where Ek c ,k t is defined in Eq. pj). For B = 0, ener- 
gies E^ kt depend only on the product rs and therefore 
coincide at both K r points. Because of the axial mag- 
netic field, this degeneracy is lifted and the spectrum 
in Eq. (11) consists of two Kramers doublets. Using 
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Eq. ( 12 1 and expanding Eq. ( 11 ) in the small parameters 
Ai/fc g < 1, k t /kg < 1, $ AB /I®oRk g ) < 1, we obtain 



E£{B) = ±hv\k g \+Ts(Ao±siga{k g }A 1 ) 
2 



1 

+ s ( -gpL B B ± ts sign{fc g }hv-^- ) , (15) 



where ± refers to the conduction and valence bands, re- 



spectively, and sign{fc s } = k g /\k g \. We will use Eq. (15) 



in conjunction with the experimental data for A$o of 
the electron and hole Kramers doublets of Fig. 4(c) of 
Ref. [HI including their ^-dependencies, in order to find 
A and Aj. This can only be done with the accuracy to 
the sign of k g . This sign remains unknown and there are 
still controversies in the literature encountering different 
theoretical models, e.g., see discussion in §5 of Ref. [9] In 
what follows, we accept k g < 0. 

Because the slopes of the ^-dependencies in Fig. 4(c) 
of Ref. [TS] are larger for the lower components of the 



Kramers doublets, both for electrons and holes, Eq. ( 15 ) 
suggests that 



ts sign(fc ff ) = ±1, 



(16) 



for the conduction and valence band, respectively. Then, 
with k g < 0, it immediately follows from Eq. (15 1 that 



2(Ao-Ai), 



a h 
^SO 



-2(Ao + Ai), (17) 



and, with the experimental values A|| = 0.37 meV 
and A§ = 0.21 meVp^O we arrive at Ai = -(A| G + 
A| G )/4 « -0.15 mcV and A = (A| G - A^ )/4 w 0.04 
meV. These parameter values will be used in all calcula- 
tions below. 

In what follows, we restrict the analysis to the lowest 
longitudinal mode. While the effect of mode mixing due 
to Coulom b repuls ion is essential in suspended wide-gap 
nanotubespSESIMl m coated narrow-gap nanotubes it is 
suppressed due to a small effective mass and large k. 



B. Nanotube Double Quantum Dots - Single 
electron regime 

In this section we outline a model of a CNT-DQD mo- 
tivated by the experimental setup in Ref. 1201 The under- 
lying geometry is sketched in Fig. [2] The potential along 
the nanotube is controlled by top gates. For simplicity, 
we choose a double-well potential to model the double 
dot, similarly to Refs. l3"4l and 1331 



V(t) = 



V , \t\<d/2, \t\ > (d + £)/2 
0, d/2 < \t\ < (d + £)/2 



(18) 



and solve the eigenvalue problem for the Hamiltonian 

H{t) = hv{Tk c a x + k t a 2 ) + ^A r + ^g^s) s + V(t), 

(19) 
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FIG. 2: Energy diagram of a symmetrical CNT-DQD with 
electrostatic gates inducing a potential V(t) along the nan- 
otube axis; the origin is chosen in the center of the double 
dot, and the energy reference point is chosen in the middle 
of the gap. Two dots D\ and D2 of length 1/2 each are tun- 
nel coupled via a middle barrier M of width d. The bound 
state energy E (shown by dashed line) obeys the criterion 
E„ < E < Vo + E g . 



with k c of Eq. (12 1. The potential V(t) is considered 
as step-like on the scale of the Fermi wavelength, 2w/kt, 
but smooth on the scale of the inverse Brillouin momen- 
tum, 2tt/\K t \. Therefore, it does not induce essential 
K + K _ -scattering. This is indeed relevant for an electro- 
statically confined dot. We calculate the single electron 
wavefunction for the symmetric geometry shown in Fig. [2] 
for equal confining potentials of the left and right wells, 
Vi = V2 = Vq. The two QDs are connected by a barrier 
of width d. The generalization to an asymmetric geome- 
try is straightforward and will result in a larger electron 
wavenumber for the deeper well. The total energy of 
the system, E = E\ + E2, will be considered fixed. In 
Eq. (19), the term in the second parentheses produces 



an (s, t) dependent level shift but does not influence the 
wave functions. For the electronic wavefunction ij) T {i) 1 
defined in different intervals, see Fig. [2] we use the ansatz 



Ae q±t ( k c-iq t 



De -1tt I Z k iq t ) +Ee 1tt ( Z k c ,- lqt 



1 J \ 1 

%U U. \ ^ —■ih,+ ( %V, 



He- qtt [ k ^ qt 



(20) 
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with z r factors defined for electrons and holes, respec- 
tively, according to Eq. Q. In the gate regions, the 
electron wavenumber is determined as 



(li 



(E ko , kt - v r/(hvy 



(21) 



Since the Hamiltonian H (t) of Eq. ( 19 1 remains invariant 
under complex conjugation accompanied by k t — > —kt, 
the wavefunctions of Eq. ( p0| can be chosen in such a 
way that C = B* , G = F* , and A, D 7 E, and H are real. 

Energy levels are found from a transcendental equa- 
tion that follows from the continuity of the wavefunction 
ip T (t) at all potential steps. Since v t = dH/dpt = V&2, 
this equation maintains the current continuity. All coef- 
ficients in Eq. (20) can be calculated from the bound- 
ary conditions and the normalization condition. For 
< Vq < 2E g , at least one electron bound state exists in 
each of the dots, and we will restrict ourselves only with 
the ground state (GS) mode. In a double dot, it splits 
into two, bonding and antibonding, with energies E\, and 
E a i,, respectively . 

The energy spectrum of the bonding mode of the 
Hamiltonian H(t) of Eq. (19 1 is plotted in Fig. [3ja) as 
a function of an axial magnetic field in the absence of 
SO coupling. Having in mind narrow gap nanotubesf^ 
we choose k g — —0.045 nm -1 , which produces a gap of 
about 2E g w 46 meV that is eight times less than es- 
timated according to Eq. ^ for v = 1 semiconducting 
tubes with R ~ 1 nm. In the absence of SO coupling, 
the B = state is four- fold degenerate in s and r, and 
its splitting by magnetic field B is shown in Fig. j3^a) . 
Fig. [3jb) demonstrates the splitting of the B = state 
by SO coupling into a doublet of two pairs of Kramers 
conjugate states, (\K+),\K^)) and (\K+),\K^)). The 
zero-field splitting is A e so ss 0.37 meV, according to the 
data of Ref. [TS] and in a reasonable agreement with the 
data of Ref. HH The B-dependence in Fig. gb) repro- 
duces correctly the patterns observed experimentally in 
Refs. EU and CS3 The effective ^-factor, dAE/d{n B B\\), 
of the upper (lower) doublet is smaller (larger) due to the 
coupling of real spin to isospin. 

The upper branch crossing in Fig. |3jb) persists for 
B || T, but turns into an avoided crossing when B ac- 
quires a perpendicular component!^ Spin relaxation near 
this crossing was investigated in Refs. l^ffllMl andflSl The 
lower branch crossing, at much lesser fields B, turns into 
an avoided crossing by K + K scattering with large mo- 
mentum transfer, A K + K -. In high quality nanotubes of 
Refs. IT81 and where both Ago and A K + K - were ob- 
served experimentally, A K + K - was nearly one order of 
magnitude less than Ago] m what follows, we disregard 
Ak+k- ■ Its origin is not well understood for now, and it 
is usually attributed to electron scattering by defectsPS 

Next, we discuss the longitudinal part of the wavefunc- 
tion ip T (t) and begin with exact symmetries of wavefunc- 
tions. Substituting Eq. (12) into Eq. (19), one arrives 
at a Hamiltonian H T:S _ B (t). It includes the parameters 
(t, s, B) only as time-inversion symmetric products tB 
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FIG. 3: (color online) Magnetic field dependence of the 
energy of the bonding mode of a single-electron DQD for 
k g = -0.045 nm" 1 and I = 200 nm, d = 120 nm, V = 3.95 
meV. The ground state energy Eos is close to the gap edge, 
Egs — Eg = 1.4 meV. (a) SO coupling is absent, Aso = 0. 
At B — 0, the spectrum is four-fold degenerate, (b) SO split 
energy spectrum. The zero field splitting between Kramers 
doublets is Aso = 0.37 meV (with Ai = —0.15 meV and 
Aq = 0.04 meV), see text for details. 



and ts. Therefore, i/_ Ti _ s ._B(^) = H TS B{t), and wave- 
functions can be chosen in such a way that they possess 
the same symmetry 



A{B Mt) 



(22) 



for an arbitrary potential V(t). The superscripts A and 
B indicate sublattice indices. For symmetric dots, V(t) — 
V(—t), one more relation holds. Performing a canonical 
transformation of i? T ,s,s(0 with a matrix <j\, one notices 
that the sign change of the term k t <J2 can be compensated 
by a t — ¥ — t transformation. Because the <Ji transforma- 

f ib A \ ( ib B 
tion transposes sublattices, o\ I ^ B J = I ^ |. wave- 
functions obey the further relations 



A{B kt) 



i ,B{A) 



(23) 



Here and in what follows we choose wavefunctions in such 
a way that they are real in the classically forbidden re- 
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gions, which is always possible because the z-factors of 
Eq. Q are real there. We note that in the left and right 
hand sides of Eq. (23) the A(B) sublattices are inter- 
changed. This reflects an absence of the microscopic 
t — > — t symmetry in chiral nanotubes. Therefore, it 
is also absent in "symmetric" DQDs despite the macro- 
scopic symmetry of the confining potential. Manifesta- 
tions of this asymmetry are discussed next. Numeri- 
cal data for the wavefunctions are presented in Fig. [4] 
Due to the exact symmetries of Eqs. (22) and (23), it is 



enough to display only a few curves demonstrating the 
basic regularities. Real parts of the wavefunctions of the 
bonding and antibonding modes are shown in Fig. (4^a) . 
Their different behavior inside the tunnel barrier is dis- 
tinctly seen. Also, A and B components of wavefunctions 
are nearly symmetric and antisymmetric for the bonding 
and antibonding modes, respectively. Technically, the 
asymmetry arises due to the admixture of the valence 
band wavefunctions and is small because the GS bind- 



ing w 1.4 meV, is small compared 
23 meV (or, what is essentially the 



ing energy, E GS 
to the gap, Eg s 
same, k tl q t <C \k g \). This asymmetry increases with 
(Ecs — Eg) and can become of the order of unity for 
\Eqs — E g \ ~ E g W^ The opposite signs of the functions 
Tp A (t) and i/> B (t) originate from k g < 0. Fig. Pub) displays 
small differences in the electron densities on Doth sublat- 
tices and their asymmetries that have the same origin as 
in Fig. [I]; a). 

Within the GS approximation, the appearance of 
bonding (lower) and antibonding (upper) tunnel com- 
ponents, denoted as ip h and ip ab , respectively, see Fig. |ij 
motivates a treatment of the orbital degrees of freedom of 
a DQD in terms of an effective two level system. Its eigen- 
states may be obtained by hybridization of the electron 
states localized in D\ or Z?2- Using the wavefunctions of 



Eq. (20) and notations of Eq. (23), we define the orbital 
basis states for certain spin and isospin (i.e., \LK T S ) and 
\RKJ)) in the left and right halfspaces, respectively, as 



\lk;) = i + < a , JJ (t)] 1 

\RKl) = i [<.,*(*) - KM*)} ' ( 24 ) 



with associated energies E TS (B); these energies include 
also the second term of Eq. (19). Each of these func- 



tions is a two-spinor in sublattice space defined by its 
A(B) components, and KJ indicates that they should 
be chosen with k c of Eq. (12) calculated for proper val- 
ues of (t, s). This parametric dependence of the orbital 
functions on s and r stems from SO coupling and will 
be essential for the classification of quantum states, see 
Sec. HIA] below. 

The connection between these basis states via tunnel- 
ing is naturally established by the overlap of the wave- 
functions in the interval \t\ < d/2. The resulting spin and 
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FIG. 4: (color online) Single particle wavefunctions i>+j? B (t) 
of a DQD for the same parameter values as in Fig. IS] and 
B = 1 T. Vertical lines sketch edges of both dots, (a) Real 
parts of i>^}^ B (t) on both A(B) sublattices are shown for 
bonding and antibonding modes, denoted by the superscripts 
a and ab, respectively. At t = the former modes retain 
a considerable amplitude, whereas the latter vanish at the 
origin. Note that at each point the A and B components 
have nearly the same absolute values but opposite signs, (b) 
Electron density distribution in the bonding state; a non- 
vanishing density at t — is visible as well as a t — > — t 
asymmetry and difference between A and B densities. 



isospin conserving Hamiltonian expressed in this basis is 



H 



RL 



Vr : s 



rprsL 

^k a ,k t 



(25) 



where E\ s \ 



\ {El 8 + E™) are the single 
particle energies and s is the detuning between the left 
and right dot. The connection between left and right 
halfspaces, established by tunneling, produces bonding 
and antibonding states of Eq. (19). The tunnel matrix 
element rj T ^ s is related to their energy difference as 

2r lT . s = El s -EH (26) 

Equivalently, the tunneling part Ht of the Hamiltonian 
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FIG. 5: (color online) Tunneling matrix elements rj TtB {B) 
plotted versus magnetic field. Their absolute values are larger 
for the lower Kramers doublet (K^ , JCjT). For B = 0, they 
depend on spin and isospin through the product sr. The zero 
field difference in the tunneling rate originates from Ai 7^ 
and vanishes for Ai — > 0. Parameter values are the same as 
in Fig.gb). 



of Eq. (25 ) can be rewritten in terms of the orbital func- 
tions \LjR)K T s ) in a form that is more convenient for 
further calculations 



H t — 



t=±,s=t,4. 



\LKT) A W(RK:\ A W + h, 



The matrix elements r\ r s are plotted in Fig. [5 
function of B for the parameter values used in Fig. 
The zero field difference in the tunneling matrix elements 
stems from the Ai term in Eq. ( 12 1. They are nearly lin- 




ear functions of B in the whole region of magnetic field 
and their absolute values are larger for the lower Kramers 
doublet. This fact is also reflected in the following in- 
equality for the wavenumbers in the classically forbidden 
regions, q t (K+) > q t (K+), cf. Eq. pB, that holds for 



our set of parameter values. Since the Hamiltonian H(t) 
of Eq. ( 19 1 is symmetric with respect to time inversion, 

r] T , s (B) = n- T ,- s {-B). 

We note that a change of the sign k g , to k g > 0, 
with the values of Ao and Ai recalculated properly to 
keep the right slopes in Fig. ^b) (A = —0.04 meV and 
Ai = —0.15 meV), keeps Fig. [5] practically unchanged 
with the upper doublet having a lesser tunneling rate 
(absolute value of 77), since the sign of Ai remains in- 
tact. Therefore, these data do not constrain the sign of 
kg either. Such a behavior of the tunneling rates seems 
counterintuitive, but it stems from the fact that the rates 
are controlled by the sign of Ai while the nature of the 
upper and lower Kramers doublets is insensitive to it. 

Fig.[6]shows the effect of detuning e on the energy spec- 
trum of the Hamiltonian Eq. (25). In our calculations, 



we have consistently used the r, s, and B dependencies 
of the matrix elements rj TiS (B) derived from Eq. (26 1. 
Cyan (dash-dotted) and red (dotted) lines in Fig. |6|a) 
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FIG. 6: (color online) Detuning dependence of the energy 
spectrum of a single-electron DQD. Tunnel matrix elements 
rj T)S from Fig. [5] are used for the respective states, other pa- 
rameters as in Fig. [3] (a) Both Kramers doublets (identified 
by color) split into bonding and antibonding modes. The 
numbers denote energy level multiplicity, (b) Level splitting 
by magnetic field B = 1 T. Because of the dominating effect of 
the Aharonov-Bohm flux, K~ components with both s =7,4- 
move upwards. The magnitudes of the avoided crossings at 
e = 0, seen both in (a) and (b), are controlled by 7j TtS . 



correspond to the upper and lower Kramers doublets. In 
the absence of a magnetic field all states are two-fold 
Kramers degenerate. For this degeneracy is lifted; 

the spectrum for B = 1 T is shown in Fig. [6jb) . En- 
ergies corresponding to various bonding (antibonding) 
and spin (isospin) states are shown in the same colors 
and line patterns as in Fig. [3j The spectrum splitting 
in a magnetic field originates from the Aharonov-Bohm 
flux &ab of Eq. ( 12 ) (that can be ascribed to an orbital 
magnetic moment [A or b) and the Zeeman spin splitting 
described by the nex t-to-last term in Eq. (19). Because 
A*orb ~ 10/xb ] 17 l 18 l 2 ° l the first contribution dominates and 



both spin components of the K~ state move up in energy 
with increasing B. The asymmetry in the level split- 
tings inside the upper and lower quadruplets in the large 
B region originates from the interference o f the SO and 
Zeeman splittings in each of the single dots) LL8122] 
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III. TWO ELECTRON REGIME 

Charge states in DQD sys tems a re usually described in 
terms of stability diagrams! 1 * 20 * 43 * Either from Coulomb 
blockade peaks in transport measurements or from charge 
sensing probes, the number of electrons and tir in the 
respective dot is monitored. In the following we consider 
the two electron regime. There are two possible physi- 
cal realizations of it depending on the adjustment of the 
confinement potentials for the left and right dot and the 
tunneling barrier between them. First, two electrons are 
confined to a single dot, here the right dot, denoted the 
(02) configuration. Otherwise, both electrons are con- 
fined in different dots, which is referred to as the (11) 
configuration. Experimentally, it is possible to drive the 
DQD between the two regimes by applied gate voltages 
if both dots are connected by a tunneling barrier. 



A. (02) Configuration 

For a two-electron system without SO coupling, a 
classification of two-electron states in terms of singlet 
and triplet states in real spin is exact due to the Pauli 
principle.^ Here, when constructing two-electron states, 
we use the Hilbert space spanned by the lowest energy 
orbitals in each dot. Because the spin and orbital degrees 
of freedom are coupled, a classification in terms of spin 
singlet and triplets is no longer applicable. Technically, 
this means that a full basis of two-particles states in the 
space of functions respecting the Pauli exclusion princi- 
ple cannot be constructed in terms of products of spin 
singlets (triplets) and linear combinations of products of 
the orbital functions \L(R)Kl) of Eq. (24 1. Spin singlets 



and triplets become coupled, and using such generalized 
basis states we arrive at the results that are quite gen- 
eral and, in particular, can be applied in the vicinity of 
the point where k c vanishes due to the cancellation of 
different terms in Eq. (12); the latter regime has been 



reported recently for a single dot However, the regime 
where the gap between the conduction and valence bands 
closes completely (see Ref. l47|) is not addressed here, since 
it does not allow for electrostatic confinement of electrons 
or holes. 

With two electrons bound to the right dot, there are 
six linearly independent antisymmetric basis functions. 
Two of them, with both spins either up or down, can be 
considered as components T±x of a spin triplet. They are 



m 2 ) = 



i 



(\RK+s} 1 \RK-s) 2 -1^2) 



(28) 



with |$" 2 ) for s =t and |$§ 2 ) for s —I, and the symbol 
1^2 indicates electron transposition. Here \RKfs) 
are products of orbital functions \RKf) and their spin 
counterparts \s). 

Four different functions, all with opposite spins, are 
spin-isospin coupled. Within the first pair of states, both 



electrons reside in the same K T point, but in such a way 
that in the B = limit one of the electrons belongs to 
the upper and the second to the lower Kramers doublet 
of Fig. [3] 



I*" 



1 
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(\RK^t)i\RKl 1)2-1^2) 



(29) 



with |$° 2 ) for r = 



and |$f ) 



for r = — . In the second 
pair of states, for B = 0, both electrons populate the 
upper and lower Kramers doublets of Fig. [3] 



\K 2 ) = \ (ml th\RK± ;> 2 - 1 ^ 2 ) . 



$ 02 ) = 



RK- t) i \RK+ i) 2 -l-2). (30) 



We designate them as 

Equations (28)-(30) demonstrate the profound effect 



|$5 2 ) and |$g ), respectively. 



of SO coupling on the symmetry of the (02) multiplets: 
Since spin and isospin are coupled, those states cannot 
be represented in terms of spin singlet and triplets. The 
energy spectrum is SO split. In the absence of a magnetic 
field, 5 = 0, the state | <& 5 ) (| $ 6 )) has the highest (low- 
est) energy because both electrons populate the upper 
(lower) Kramers doublet, while the four other states are 
mutually degenerate because one of the electrons popu- 
lates the upper and the second the lower state. 



B. Coulomb Matrix Elements 



We choose the Coulomb potential as 



Vn 



yfit + (h - t 2 ) 2 + {2R? sin 2 ^ - ip 2 )/2] ' 

(31) 

with k as an effective dielectric constant. The cut-off 
term a 2 (2as) 2 , with clb for the Bohr radius, accounts 
for the size of 2p z functions.^ Such a cut-off is convenient 
in numerical calculations but has no essential effect on 
the final results because Coulomb matrix elements con- 
verge in two dimensions at small electron separations. In 
what follows, we calculate matrix elements of Vc for both 
the (02) and (11) configurations. Similarly, while taking 
into account consistently the SO corrections to functions 
\KJ) in both dots, we keep only spin-diagonal terms when 
calculating the Coulomb matrix elements and therefore 
exclude spin nonconserving processes. This approxima- 
tion is motivated by our focus on dots with small radius 
R and narrow gap E g w hvk c ; indeed, the diagonal SO 
corrections are large in inverse k c while nondiagonal SO 
terms are suppressed for small R by strong orbital quan- 
tization in the circular directionPZl In exchange matrix 
elements for the (02) configuration, a selection rule for r 
appears from the fact that Vc depends on (^1,(^2) om y 
through their difference. Upon using notations |rs) for 
kets including products of the orbital functions \KJ) and 



10 



the corresponding spin functions |s), we take advantage 



of Eq. ( 13 I and find 



(riSi(l),75« 2 (2)|V r c (l,2)|7S«3(2) ) r 4 S4(l)> 

/-2tt 

cx / dipexp [—i(M + n 2 /2)(ri + t 2 - r 3 - t 4 )c^], 
Jo 

(32) 

where if = (tpi + ip 2 )/2. Expressing M in terms of the 
chiral indices (ni,n2) results in (M + n 2 /2) = m/3 + 
n 2 /6 — i^/3. Hence, inside the irreducible wedge of the 
Bravais lattice where < 6 < tt/6 and ni,ri2 > 0, it is 
always true that (M + n 2 /2) ^=0, and therefore 




Tl+T-2 = T 3 



r 4 . 



(33) 



s 4 , 



Comparing Eq. (33 1 with the spin selection rule Si 
S2 = S3 (hi the leading approximation in SO) underscores 
a fundamental difference between the spin and isospin, 
which is an orbital quantum number. 

Under these assumptions, Vc of Eq . |31~1 ) is represented 
in the basis of <& 02 -functions of Eqs. (|28|)-(30) as 



J, 



t 



112 



V 



u, 



(12 
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(34) 

It includes six Coulomb matrix elements Uq 2 and three 
exchange matrix elements Jo 2 . The latter ones include 
terms that are nondiagonal in r but obey the selection 
rule of Eq. (33); one can show that the nondiagonal ma- 
trix elements Jo 2 are real. The absence of isospin indices 
in Uq 2 and Jq 2 indicates that electrons belong to differ- 
ent valleys. Similarly, the absence of spin indices in Uq 2 
indicates that electrons possess opposite spins. 

All Coulomb terms U02 have a universal form in the 
framework of our model and do not depend on the chi- 
rality of the nanotube. As distinct from them, exchange 
integrals Jo 2 depend on chirality and, what is even more 
important, include products of orbital functions \RKJ) 
with different values of r. Therefore, they require large 
momentum transfer of 47r/3a, which is accompanied by 
fast oscillating factors in the integrands. The calcula- 
tion of such integrals cannot be performed using envelope 
functions of Eq. ^ and requires including microscopic 
Bloch functions of graphene and short range interaction 
potentials. This is outside the framework of our model, 
and since such integrals are small, we disregard them in 
what follows. 

The six different Coulomb integrals Uq 2 and their mag- 
netic field dependence are shown in Fig. [7] for k = 1. As 
distinct from the single electron levels of Fig. [3j where 
the B dependence originated from the orbital and spin 
magnetic moments, the B dependence of Uq 2 integrals is 
determined by the B dependence of the circumferencial 



FIG. 7: (color online) Magnetic field dependence of the 
Coulomb matrix elements U02 in a single QD with dielectric 
constant k — 1. Parameters are the same as in Fig.[3fb). The 
B dependence is strong for isospin polarized statesT^ 2 ) and 
|$2 2 ), with the opposite signs of the slope for r = ±. For 
isospin unpolarized states the B-dependences are weak. The 
energies Uq 2 an d Uq 2 nearly coincide. For B = 0, the largest 
difference in U02 energies is achieved for the states with both 
electrons belonging to the lower or upper Kramers doublet. 
The matrix elements C/02 strongly influence the position of 
the (ll)-(02) degeneracy point in Figs.[9]and 10 below. 



wavenumber k c of Eq. ( 12 ) only. Remarkably, this de- 
pendence is much stronger for isospin polarized states 
I $§ 2 ) and I $4 2 ) than for the other four states; for 
the latter ones, the B-dependencies are nearly identi- 
cal. We attribute this behavior to the competition of 
the two largest terms in fc c , namely the first and sec- 
ond one of Eq. (12), which therefore does not rely on 
SO coupling. Iny~3>3 2 ) and | <I>4 2 ) both electrons have 
the same isospin t, hence the same B-dependences of k c 
add, while in all other functions the electrons have oppo- 
site signs of r and the Independences subtract. We note 
that in the absence of SO coupling, Coulomb integrals 
for I $ ( j )2 ), I c & 2 2 ), I <1>5 2 ), I ®q 2 ) coincide for all magnetic 
fields. For nanotubes coatedby an insulator, as in ex- 
periments by Churchill et alP23, the Coulomb repulsion 
is reduced by a factor of n ~ 10. It should be noted that 
the metallic gates used in the experimental setups also 
strongly screen the Coulomb interaction, and in particu- 
lar they cut-off the long range part of it. Therefore, the 
absolute numbers for the Coulomb matrix clement that 
we find here are subject to changes depending on the ex- 
perimental details, however the general trends based on 
the symmetry of the two-particle wavefunctions remains. 



C. Energy spectrum 

The results for the magnetic field dependence of the 
energy levels of a two-electron QD are shown in Fig. [8] 
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FIG. 8: (color online) Magnetic field dependence of the en- 
ergy spectrum of a two-electron single QD. Coulomb interac- 
tion is screened by k = 10; other parameters are the same as 
in Fig. [3] (a) Without SO coupling the spectrum is six-fold 
degenerate at B — 0, and its B dependence originates mostly 
from the coupling of the orbital and spin magnetic moments to 
the field. Wavefunctions can be represented as spin singlets- 
isospin triplets and spin triplets-isospin singlets, (b) With SO 
coupling the spectrum is split at B = 0. The level crossing at 
finite B results in a ground state change from two electrons 
populating the lower Kramers doublet to two isospin polar- 
ized electrons. Numbers with arrows denote the energy that 
corresponds to a particular state among | ^if...^)- 



We have diagonalized the two particle Hamiltonian for 
a single QD, cf. Ref. [331 ln the presence of SO coupling 
as well as a screened Coulomb interaction, k — 10. QD 
parameters are chosen as in Sec. |IIB| The B dependence 
of the Coulomb interaction terms Uq2 was taken into ac- 
count consistently. The comparison of panels (a) and 
(b) demonstrates the effect of SO coupling. In the ab- 
sence of SO coupling, Fig. [8^ a), the hexaplet remains 
unsplit at B = because J02 exchange integrals are dis- 
regarded. For A50 = 0, a classification of these degener- 
ate states in terms of a spin singlet (isospin triplet) and 
a spin triplet (isospin singlet) is appropriate. The B de- 
pendence of the isospin polarized states l^ 2 ), l^ 2 ) and 
spin polarized states l^? 2 ), l^" 2 ) originates mostly from 



their orbital and spin magnetic moments, respectively. In 
the whole region of magnetic fields, the ground state is 
isospin polarized with both electrons in the t — 1 state 
having opposite spins. SO coupling splits the hexaplet 
at B = 0, see Fig. |8^b), with both electrons populat- 
ing the lower Kramers doublet in the lowest state I'&ep)- 
The level crossing at B sa 0.3 T results in a change of 
the ground state. This transition might also be seen in 
the two electron spectrum of Ref. HH1 At larger fields, 
the ground state becomes well separated from all higher 
states. 

Although we do not calculate electron attraction due 
to their coupling to phonons, we note that an estimate 
shows that it might become comparable to a screened 
Coulomb repulsion for k > 10. A more detailed investi- 
gation of this contribution is needed. 



IV. TWO PARTICLE SPECTRUM AS A 
FUNCTION OF DETUNING 

In recent experiments,^ the dephasing time T 2 * of a 
two particle state was obtained by the following mea- 
surement cycle. First the system is prepared in the (02) 
configuration whose ground state is non-degenerate for 
a CNT-DQD, Fig.^b). The doubly occupied right dot 
might be considered as a double dot in a strongly de- 
tuned state where the detuning energy e compensates the 
strong Coulomb repulsion; hence, it is energetically favor- 
able for two electrons to populate the same dot. When 
decreasing e, the Coulomb repulsion and interdot tun- 
neling allow pushing one of the electrons to the left dot, 
and the system is transferred into the (11) configuration. 
This produces an additional degree of freedom, manifest- 
ing itself in a quantum well index L(R), and allowing 
for 16 states in the (11) configuration, as compared to 
6 states in the (02) configuration. The whole space in- 
cludes 22 basis states. The transfer of the system from 
the six-fold (02) space to sixteen-fold (11) space is fol- 
lowed by dephasing due to different mechanisms, includ- 
ing hyperfine interactions and isospin scattering. When 
e is increased again, after a certain separation time T5, 
the system is prevented from coming back because not all 
states from (11) are connected by tunneling to the states 
in (02). This generalized Pauli blockade arises from the 
selection rules both in spin and isospin. The probability 
of finding both electrons again in the right dot depends 
on rs, and measuring the return probability as a function 
of ts is used to extract T, , as has been done for CNT- 
DQDs with T* = 3.2 nsPf Therefore, T* should strongly 
depend on the coupling between the energy levels of the 
(02) and (11) subsystems that will be investigated below. 

A. (11) configuration 

In this section we construct a basis of the two-particle 
Hilbert space of the (11) configuration starting from 
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the (02) configuration basis and using the correspond- 
ing single-particle wavefunctions of Eq. (24). Consider, 
e.g., the state of Eq. (28) with spin polarized functions 
| Q® 2 )- There are two possibilities, either the first or the 
second electron can occupy the right dot, 



i 



(| RK+a)i I LK-s) 2 - 1 ^ 2) , 

(| LK+s)i | RK-s) 2 -1^2). (35) 



This procedure, when applied to the states 
Eqs. (p8l, p9l) and pOl 



(02) 



of 



1,...,6 

results in twelve states of the 



(11) configuration. From those twelve states we construct 
combinations that are symmetric and antisymmetric in 
L/R space. For example, from Eq. (35) we find for s =f 



1 
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± 



if 



(36) 



We denote by | <5>\\ 6+ ) and | 



,...,6- 



respectively, 

the symmetric and antisymmetric superpositions in L/R 
space. In addition, in the (11) configuration there are 
four states polarized in both spin and isospin 



1 



LKIf 



RKls) 



1 



(37) 



All of them are antisymmetric; similar combinations in 
(02) are forbidden by the Pauli exclusion principle. We 
use the convention i = 13 for t — +,s i = 14 for 
t = +,s =^, i — 15 for r = —,s =t, and i = 16 for 

T = -,S =1. 



B. Coulomb and tunneling matrix elements 

We need to calculate a 16 x 16 matrix of the Coulomb 
interaction in the (11) configuration that is similar to 
Eq. (34), as well as one- and two-particle matrix elements 



that connect the (02)- and (ll)-subspaces. 

In the (11) subspace, the matrix of direct Coulomb 
terms is diagonal and its matrix elements are equal 
for symmetric and antisymmetric combinations. Hence, 
we need to compute only six independent matrix 
elements for the twelve symmetric and antisym- 
metric states |$i± 6±)- They are denoted as 
U^ 1 ,U n , 11^,11^,11^,11^ according to the spins and 
isospins of the states involved. There are four additional 
Coulomb terms for the states | 14 15 16 ) that are both 
spin and isospin polarized. We denote their Coulomb 
matrix elements as lf[± , U^, lf\\ ■> • Since we have 
chosen the wavefunctions in such a way that they are 
real in the classically forbidden regions (see Sec. II B 
Coulomb matrix elements between the states | <&i +j ... ! 6+, 
and | $1—^ e— ) vanish. Furthermore, Coulomb integrals 
between |3>i±,...,6±) and | < E , i3,i4,i5,i6) vanish because of 
the spin and isospin selection rules. 



The strongest ^-dependence occurs for the six matrix 
elements ( U^ x , U xl , , , Jj\^ , llf^ ) corresponding 
to the isospin polarized states (not shown here). Simi- 
larly to the (02) configuration, we attribute this behavior 
to the -B-dependence of k c of Eq. ( fl2| ). However, compar- 
ing to the (02) configuration (see Fig. [T| the slopes for 
and | $3±) have opposite signs. The same holds 
and | 

In addition to the matrix of the direct Coulomb inter- 
action discussed above, there are also exchange matrix 
elements. As distinct from the (02) configuration, in the 
(11) configuration there exist a number of interdot ex- 
change matrix elements that are not annihilated by the 
requirements of the spin conservation and the selection 
rule of Eq. ( 33 ) . They include the overlapping densities 

U \ 

!+/> 



also for | 



between the left and right dot. Specifically, for | &} 1 



Jl R = (LK-2,RK. 



t,l 



V c 



LK+^RK-J 



(38) 



The generalization for other states is straightforward. We 
find that there are 10 independent J-exchange matrix ele- 
ments in the (11) configuration. All of them arc positive, 
and they are one order of magnitude smaller than the 
Coulomb terms for the parameters values chosen. In the 
(11) configuration, these exchange terms shift the ener- 
gies of the antisymmetric states down and the symmetric 
states up. 

Besides the (11) terms in the Hamiltonian matrix cal- 
culated above, the cross-terms that provide a coupling 
between the (11) and (02) configurations are of critical 
importance. They originate both from the single-electron 
tunneling Hamiltonian Ht and from the two-electron 
Coulomb Hamiltonian Vc and connect all |<I> 02 ) states 
with the first twelve l^ 11 ) states. The last four |$ n ) 



states of Eq. (37) cannot tunnel to the (02) configuration 
by construction. 



The first contribution is similar to Eq. (25), yielding 
twelve matrix elements 



$° 2 I H T I $ 



s±l 



11 \ 

T±) 



1 

i 



{V+,s ±77_, 



($" 2 I H T I $ n 



u±l 



($° 2 I H T I 



V2 



(»7-4 ± »7+,t) . 
0?+,4- ± V-,t) ■ 



(39) 



Here Ht are sums of the tunnel Hamiltonians of Eq. ( 27 ) 
over both electrons, and rj T S are defined by Eq. (26). 



Note that ± signs in Eq. ( 39 ) correspond to symmetric 
and antisymmetric states in the (11) configuration. At 
B = 0, antisymmetric combinations |5?5_ g_) are forbid- 
den from tunneling to (02) because of spin and isospin 
conservation and the relation rj T S { B) = r]— T ,— g (—B) es- 



tablished in Sec. II B From Eq. ( 39 1 follows that antisym- 
metric combinations are not entirely forbidden from tun- 
neling to the (02) configuration because of the (r, s) de- 
pendence of the tunneling matrix elements r/ r s . However, 
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as one concludes from Fig. [5] this dependence is rather 
weak, only about 1%, and therefore transitions to (02) 
states from antisymmetric (11) states are strongly sup- 
pressed. It is in this sense that we discuss the left-right 
(L/R) symmetry selection rules in what follows. Since 
they are approximate, the corresponding level crossings 
transform into narrow avoided crossings. 

The second contribution originates from the Coulomb 
interaction and is also represented by 12 matrix elements 



<d>° 2 | V C | 

(K 2 1 v c I K±) 



7 s 
j Q 



(0,02 



v c I = Ji 



s± 

02-f+lU 

T± 
TU± 

J 02-H-ll' 
l± 

02-f+ll- 



(40) 



All of them include overlap densities from the right and 
left dots. As an example, we present the exchange inte- 
gral between the states |$? 2 ) and |$}±) 



4Ln = -j= {(RK^RK+i | V C | RK+iLK^) 

± (RK^ 2 RK+ a | Vc | LK+^RK^)) . (41) 



Other exchange integrals of Eq. ( 40 1 have a similar struc 



ture and the ± signs refer to the symmetric and antisym- 
metric linear combinations in the L/R degree of freedom 
for the (11) configuration. The integrals of Eqs. (39) 



and (41| are subject to the same spin/isospin selection 



rules and contribute additively to all (avoided) crossings 
between the (02) and (11) states. Comparing Eqs. (38) 



and (41) one notices immediately that matrix elements 



including an odd (even) number of the wave functions of 
the left or the right dot have opposite (same) signs in the 
symmetric and antisymmetric states. 



C. Energy spectrum 

For electron spin dynamics as well as for electron ma- 
nipulation by gate voltages, the dependence of the en- 
ergy levels on the detuning e between the two dots and 
the magnetic field B is important. Especially, the po- 
sition of the energy levels and the width of the avoided 
crossings that appear due to tunneling and exchange in- 
tegrals might be observed in transport experiments on 
CNT-DQDs. 

Fig. [9] presents the result of the two particle spectrum 
as a function of detuning at B = 0. Other dot parameters 
are the same as in Fig.[3jb). In the lower right corner, the 
system is in the (02) configuration and the ground state 
is given by | <I>g 2 ) with both electrons populating the 
lower Kramers doublet. The next group of states (l*^ 2 ) 
with s =t,4- and \§® 2 ) with r = ±) originates from the 
mixed populations of the two Kramers doublets. The 
splittings between their energies are controlled by the 
matrix elements Uq 2 and Uq 2 , the first of which is small 



and the second vanishes at B = 0, see Fig. [7J They are 
not resolved in Fig. [9| 

When e decreases, the six (02) states hybridize with 
their counterparts and form lower and upper 

(bonding and antibonding) tunnel components, indicated 
by B and AB superscripts in the figure. The (11) — (02) 
degeneracy point is located at s w U02/2. Here U02 is a 
mean value of the integrals Uq 2 with i = (f, 4-, +, — , u, I) 



defined in Sec. IIIB that differ only within 10% among 
each other. More accurate positions of the degeneracy 
points for specific transitions are e° = [Uq 2 — (U^ + 
J l LR )}/2 where J l LR are exchange integrals defined by 
Eq. (38); J\ R are only about a few percents of Uo2- The 
widths of the tunnel doublets at e° points can be found 
from Eqs. (f39j) and ((41]), for instance, for the up-spin 

As seen 



state |$i) it equals \/2(?7+t + V—f) + ^^i. 
in Fig. [9j the splitting between bonding and antibonding 
states becomes large in the vicinity of the (ll)-(02) de- 
generacy point and competes with the SO induced split- 
ting. 

We note that these equations do not involve the states 
I ^13 14.15 16) since they are completely decoupled from 
all different states and pass through the whole region of 
the (02) — (11) resonances without any avoided crossings. 

Remarkable properties of the ground state deserve a 
special attention. In the absence of Coulomb interac- 
tion (or for very large k), the ground state of the (11) 
configuration is controlled by tunneling and is a bonding 
state that is always symmetrical, see Fig.|4^a). However, 
because of the competition between tunneling, Eq. (39), 
and exchange, Eqs. (38) and (40), the ordering of lev- 
els can change. This reordering of levels is a real con- 
sequence of our calculations despite the large value of 
k = 10. The level of the symmetrical hybridized state of 
|$gi) and l^g 12 ) (designated as a bonding state in 
Fig.pl crosses the level of the antisymmetric state |$gi) 
at the point highlighted by a circle in Fig. [9j While to 
the left from the circle these levels nearly merge in Fig. |9j 
they are well resolved in Fig. 11 Under such conditions, 



the ground states on the left and right from the (11)- 
(02) resonance are not connected because of the highly 
unusual order in which the levels follow on the left, i.e., 
in the mostly (11) configuration. There the antisymmet- 
ric state lies below the symmetric one, as a consequence 



of the fact that the exchange integrals (see Sec. IV B) 
prevail over the competing contribution of the tunneling 
matrix elements. 

Besides the state |4>g-) there are different (nearly) un- 
connected states showing up as lines monotonously in- 
creasing with e in Fig [9} Altogether, there are 10 states 
in (11) from which electrons cannot tunnel to (02), six 
antisymmetric states and four states |$}| 14 15 16 ) which 
have no counterparts in the (02) configuration. 

In Fig. |9j there are two kinds of level crossings. All 
crossings related to I < &i3 j i4 i i5 ; i6) are robust, in the frame- 
work of our scheme, against small perturbations because 
these are the only states that are both spin and isospin 
polarized. In particular, they do not rely on the L/R 
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symmetry. Distinct from them, the crossings involving 
1^5 6-) states and narrow anticrossings involving other 
states rely on the L/R symmetry (that is not exact 
and is based on the weak dependence of r) T S on (r, s) and 
similar properties of exchange integrals, cf. Sec. |IV| B). A 
violation of this symmetry transforms them into avoided 
crossings, therefore, the widths of the anticrossings can 
be controlled by the gates. This might allow control of 
the system passage across the point indicated by a circle 
in Fig. [9] when performing the (02) to (11) excursions. 
Note, the above statements relate only to the stability of 
the crossings. The very fact of the existence of specific 
crossings depends on the relative magnitude of a number 
of different Coulomb and tunnel matrix elements. 
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FIG. 9: (color online) Two particle spectrum at B = as 
a function of detuning e demonstrating a gradual transition 
between the (11) and (02) configurations. Parameter values 
are the same as in Fig. |3|b), dielectric constant k — 10. Hy- 
bridized bonding and antibonding states are designated as 
$f = a i $}l+p i $'} 3 and $f s = a^Jj - ft$f , respectively. 
The coefficients cn,/3i > depend on e and were found from 
numerical diagonalization. Antisymmetric |$ lx ) states that 
practically do not hybridize with |<1? 02 } states are designated 
as I Slashes indicate the states that are either exactly 
(Kramers) or nearly degenerate; all of them are spin (isospin) 
polarized. The states | <C»5 ) and are split by SO coupling. 
Remarkably, the ground states of (02) and (11) are not tunnel 
coupled, and the circle highlights the crossing point between 
the states |$f ) and | <E>qL > 



One general comment regarding the ground states of 
two-electron D QDs i s relevant. According to the Lieb- 
Mattis theorem ] 49 ! 50 ! the ground state of a two-electron 
system, at B — and in the absence of SO coupling, is 
always a spin singlet. The proof of this statement (at- 
tributed to Wigner in Ref. 111?]) is applicable only to scalar 
wave functions. Therefore, it is not applicable to carbon 
nanotubes where the wave functions are spinors in the 
pseudospin space. Also, the classification of the quantum 
states of SO coupled systems is generically impossible in 
terms of the spin eigenstates. Hence, it is quite remark- 



able that despite all these odds, both GS wave functions 
of Fig. 9 belong to the group of functions with zero mean 
value of the spin. A specific GS function is chosen by a 
number of competing parameters. 

Fig. [lO] presents the energy spectrum as a function of 
detuning e for a magnetic field B = 1 T. All degenera- 
cies, both in (11) and (02), are lifted. This field is large 
enough to change the symmetry of the ground states both 
in the (11)- and (02)- configurations. Once again, a GS 
to GS transition is not allowed. At the (02) side, the 
splitting of the isospin doublet | < 1 >0 - 2 ), t = ±1, becomes 
larger than the SO splitting separating \&® 2 ) and |$g 2 ) 
states and its r = + component |$3 2 ) shifts to the spec- 
trum bottom, in agreement with the experimental find- 
ings of Refs. [T51 and [2TJI At the (11) side, the magnetic 
field splits the \<f>{\ 15 ) Kramers doublet (that was only 
slightly above the ground state in Fig. [9]) and shifts its 
I $14) component to the spectrum bottom; it is spin and 
isospin polarized with s =4,, r = +. In this context, it is 
instructive to follow the adiabatic evolution of the |$f ) 
ground state starting from \Q® 2 ) in the lower right corner 
of Fig. [ToJ After crossing the ground state (this 

crossing is both spin conservation and L/R symmetry 
protected), it passes through a narrow anticrossing with 
1^3-) (protected by the weak B dependence of r] TyS (B) 
and highlighted by a circ le) to appear only slightly above 
see Fig. 



it as 

ment to Fig. [9] the 
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Similarly to the related com- 
cTth of the avoided crossing can be 
enhanced by producing asymmetry between the left and 
right dots. Since \$\\) and both states |$3±) possess the 
same pseudospin r = + while \&\±) are spin unpolarized 
and 1$ 



14 / is spin polarized, the relaxation from |3>3±) 



to 



the ground state is only possible due to the spin noncon- 
servation. Therefore, excursions from (02) to (11) can be 
used for measuring the spin relaxation rate. 

We have checked the behavior of the level crossings, 
highlighted by circles in Figs. |9"[and[T0} when we change 
the size of the gap E g . With E g increasing twice, both 
crossings remain stable and move to the right, nearly half 
way to the (ll)-(02) degeneracy point (e ~ 11 meV). 

In Figs. [9] and [lOj in the vicinity of the (ll)-(02) de- 
generacy point, gross features are dominated by the \^ n ) 
- |$ 02 ) hybridization. To illuminate different properties 
of the spectrum, its SO-coupling and independence, we 
present in Fig. [TT] the energy spectrum at e — 2 meV. 
While it was found by the same procedure as Figs. [9] and 
[TO} we checked that it is very close to the spectrum found 
in the 16 x 16 basis of |$ 1:L ) functions; in particular, all 
levels follow in the same order. This proves that the 
contribution of the polar configuration (20) (with both 
electrons on the left dot) not included in our calculations 
is small at e = 2 meV and has only minor effect on the 
results. 

At B = 0, the spectrum is dominated by the split- 
ting originating from the one- or two-fold population 
of the upper and lower Kramers doublets separated by 
A50 ~ 0.4 meV. Splittings from the inter-dot exchange 



matrix elements are lesser: J 



LR 



0.03 meV for k = 10, 



15 
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FIG. 10: (color online) Same as in Fig.[9]for a magnetic field 
B = 1 T. Again, bonding and antibonding states are denoted 
as <3>f and <&f B , respectively. All degeneracies are removed. 
Four <E>i3_ ie ) states and six that practically do not 

hybridize with |$ 02 } states are shown as ascending lines. The 
circle highlights the intersection of the |$^} hybridized state 
with \<&l-)- See text for details. 



Fig. [IT] The tunneling matrix elements r\ ~ 0.02 meV 
[see Fig. [5] and Eq. ( 39 )] also induce lesser splittings. 
Therefore, the gross structure of the energy spectrum 
is controlled by A50, and this suggests describing it pri- 
marily in terms of Kramers doublets rather than inde- 
pendent spin and isospin populations. The fine struc- 
ture inside each group (4+8+4), originating from tun- 
neling and Coulomb terms, should be accessible for ex- 
perimental resolution. One can also distinguish energy 
differences between the bonding and antibonding | $,-_J_) 
states and the states forbidden for tunneling to (02). It is 
seen that the energies of antisymmetric states are lower 
than the energies of the corresponding symmetric states 
in the whole range of magnetic fields. This is the result of 
the Coulomb interaction that favors antisymmetric states 
prevailing over tunneling that favors symmetric states. 

The ^-dependence is dominated by the isospin Zeeman 
coupling because fi or b hb- However, more careful ex- 
amination allows distinguishing differences in the slopes 
of the states with the spin and isospin polarized in the 
same or in opposite directions, e.g., \&\\ ,«). A Zeeman 



splitting of spin polarized states \$\± 2 ±) is distinctly 
seen. 



V. SUMMARY AND DISCUSSION 

We have studied the detailed structure of the energy 
spectrum of a symmetric carbon nanotube double quan- 
tum dot with either one or two electrons confined by an 
electrostatic potential. We focused on narrow-gap coated 
nanotubes allowing efficient gate control for electronic 
and spintronic applications, and investigated the effect 



FIG. 11: (color online) Magnetic field dependence of the en- 
ergy spectrum of a two-electron double dot for e — 2 meV 
where the admixture of both (02) and (20) configurations 
is negligibly small. The Coulomb interaction is screened by 
k — 10 and other parameters are the same as in Fig. [3] Num- 
bers with arrows denote particular | 3 > i± l ... l 6±,i4,....i6) states. 
Relative magnitudes of the different level splittings originat- 
ing from the tunneling and Coulomb interaction are distinctly 
seen. At B — 0, the dominating splitting comes from Aso 
depending on whether one or both electrons belong to the up- 
per or lower Kramers doublet; other contributions are smaller 
by one order of magnitude. The strong B dependence is con- 
trolled by the isospin through /i or (, and a weaker one by the 
spin through {1b ■ The splittings of bonding and antibonding 
levels are weak, and their sign is controlled by the prevalence 
of the Coulomb contribution over tunneling. 



of both SO coupling constants, Ao and Ai, on the en- 
ergy spectrum. The large effective dielectric constant of 
such nanotubes (k ~ 10) in conjunction with a small elec- 
tron effective mass m* ~ E g /v 2 suppresses admixture of 
the higher longitudinal modes and allows studying the 
fine structure of the spectrum originating from the spin 
and isospin degrees of freedom in the framework of a 
single-mode theory. The importance o f suc h a study is 
called for by the experimental discover y^ 20 ! of very nar- 
row features (< 10 mT) in the magnetotransport spectra 
of DQDs. While a recent theory^ proposed a mecha- 
nism for developing magnetocurrent minima with a width 
Aso/Horb ~ 100 mT based on the global width Aso of 
the SO split spectrum, unveiling the nature of the nar- 
row features seem to require mechanisms involving spe- 
cific quantum levels. Note that the basic elements of our 
theory are also applicable to suspended semiconducting 
nanotubes as well, but accounting for higher longitudinal 
modes might become necessary! 35 ! 36 * 44 ! 

After solving a spinor equation for a double square- 
well confining potential in the axial direction, we ob- 
tained the single-particle spectrum in the presence of 
SO interaction. Due to the coupling between spin and 
isospin, the four-fold degeneracy is lifted at zero mag- 
netic field, B — 0, which results in two Kramers doublets 
corresponding to either aligned or anti-aligned spin and 
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isospin. We note that while the diagonal Ao and nondi- 
agonal Ai SO coupling constants combine in the splitting 
A50 between the Kramers doublets, Ai contributes in- 
dependently to the interdot tunneling rate. As a result, 
Kramers doublets acquire different B-dependent tunnel- 
ing rates, Fig. [5j we estimated this difference using the 
realistic values of A and Ax found from the experimen- 
tal data of Ref. 15] They can be observed in experiments 
on single-electron transport across double dots. 

The basis states for a two-electron DQD in the (02) 
configuration (both electrons on the right dot) include 
two spin polarized as well as two isospin polarized func- 
tions, and two functions belonging to the upper and lower 
Kramers doublet, respectively. All of them are spin- 
isospin coupled, and the Coulomb interaction energies 
depend both on the spin and isospin. In the (11) con- 
figuration of a symmetric DQD, these functions generate 
(by moving a single electron from the right to the left 
dot) 12 basis functions of which 6 are symmetric (bond- 
ing) and 6 antisymmetric (antibonding) in the indices of 
the left (L) and right (R) dots. Four more states, all L/R 
antisymmetric, have no analogs in the (02) space. Only 
bonding modes strongly hybridize with (02) states. 

Our main result is the energy spectrum of a two- 
electron DQD, calculated for the regime of comparable 
tunneling and SO energies, shown in Fig.|9]for B = and 
in Fig. [TO] for B = 1 T as a function of the detuning e 



between left and right dots. It is discussed in Sec. IV C 



Figures [9] and [To] illustrate how fundamentally the isospin 
degree of freedom and its coupling to the spin change the 
spectrum. This change makes the analysis of the spec- 
trum much more complicated compared to the spectrum 
of GaAs DQDs^ which consists of the spin singlet and 
triplet branches alone. 

While both the Pauli blockade and dephasing rate are 
challenging goals for experimental studies, investigating 
the dephasing rate by initializing the system in the (02) 
configuration and making excursions into the (11) config- 
uration is more tractable from a theoretical point of view 
because of a lesser manifold of quantum states whose 
width can be controlled by gate potentials. The effect of 
a magnetic field on the mutual position of the lower lev- 
els that influences the relaxation rate between them can 
be inferred from Figs. [9] and 10 as discussed in Sec. IV C 



As distinct from GaAs where the ground state is a sin- 
glet, in nanotube DQDs this is a double-populated lower 
Kramers doublet. We have also found that in our pa- 
rameter range the ground state in the (11) configuration 
is antisymmetric in L/R, indices because the Coulomb 
repulsion prevails over tunneling. This unique situation 
results in the opposite L/R parity of the ground state 
on both sides of the (ll)-(02) degeneracy point, Fig. [9] 
Therefore, low energy excursions into the (11) config- 



uration can probe the relaxation rate at small energy 
transfers and indicate the position of the L/R symmetry 
point (deviation from it turns the level crossing into an 
anticrossing). When the magnetic field becomes strong 
enough, Zeeman splitting shifts a spin polarized state to 
the bottom of the spectrum. As a result, ground states 
on the left and right differ not only in the L/R symmetry 
but also in the two-particle spin wavefunction, Fig. |10| 
Hence, similar excursions can probe the spin relaxation 
rate t s _1 . Moving up in energy should allow probing 
higher states of the (11) configuration. 

With such a rich energy spectrum, the very notion 
of the spin (Pauli) blockade should be generalized 
including both spin and isospin, and the blockade be- 
comes rather sensitive to the parameters of the system. 
Therefore, it is natural that the blockade has either been 
observed^ or alternatively not observed^ by different 
experimental groups. The outcome should strongly de- 
pend on populating the different (11) levels during the 
initiation phase, mechanisms of the relaxation and leak- 
age, and the topology of the dense intertwined net of the 
energy levels. A significant challenge is establishing the 
optimal conditions for achieving the Pauli blockade. 

The pattern of the energy spectrum, which is rather in- 
volved even in the framework of a simple model (Fig. 10 ) 
should become even more complicated in realistic sys- 
tems due to the r non-conservation that is usually con- 
trolled by extrinsic mechanisms and, therefore, might be 
different in the left and right dots. It can be taken into 
account either phenomcnologi cally by in cluding a term 
Akk'Ti into the Hamiltonian) 20 | 34 | 4 -Sl^l or by modeling 
a short-range disorderPsl Likewise, the electron attrac- 
tion through their coupling to stretching phonons that 
can compete with the Coulomb repulsion at n ~ 10, cf. 
Sec. IIII CI is not studied here and deserves a detailed 



investigation in the future. 

Note added: While completing this manuscript, we 
became aware of a paper by v. Stecher et alP^on a related 
subject. Both studies are complementary, since Ref. 0U 
focuses mostly on the effect of electronic correlations in 
suspended wide-gap nanotubes, while we concentrate on 
the fine SO structure of the spectra of coated narrow-gap 
nanotubes where such correlations are suppressed. 
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